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THE GEODESIC FLOW ON NILMANIFOLDS 


ALEJANDRO KOCSARD, GABRIELA P. OVANDO, AND SILVIO REGGIANI 

Abstract. In this paper we study the geodesic flow on nilmanifolds equipped 
with a left-invariant metric. We write the underlying definitions and find general 
formulas for the Poisson involution. As an example we develop the Heisenberg Lie 
group equipped with its canonical metric. We prove that a family of first integrals 
giving the complete integrability can be read off at the Lie algebra of the isometry 
group. We also explain the complete integrability on compact quotients and for 
any invariant metric. 


1. Introduction 

Given a smooth manifold M, any complete Riemannian structure (•, •) induces 
the geodesic flow T: M x M —>■ M which can be defined as the Hamiltonian flow 
associated to the energy function E{v) := l/2{v,v) on TM. 

Usually, this flow is not integrable in the sense of Liouville and it is generally 
expected that the integrability of the geodesic flow imposes important obstructions 
to the topology of the supporting manifold. 

However, contrasting some results of Taimanov [I9l[20] on topological obstruc¬ 
tions for real-analytic manifolds supporting real-analytic integrable geodesic flows 
with some smooth examples of smoothly integrable geodesic flows on manifolds that 
do not satisfy the above obstructions constructed by Butler [3], and Bolsinov and 
Taimanov [2], we observe that the regularity of first integrals plays a fundamental 
role that nowadays is not completely well understood. 

For that reason, when dealing with locally homogeneous manifolds to have the 
possibility of constructing real-analytic first integrals is very desirable. 

In the advances reached in the theory of Hamiltonian systems in the 1980’s one 
can recognize the role of Lie theory in the study of several examples. This is the 
case of the so known Adler-Kostant-Symes [iiiiaiiB] scheme used for the study of 
some mechanical systems and of the so known Thimm’s method for the study of 
the geodesic flow m- In both cases the main examples arise from semisimple Lie 
groups. These results appeared parallel to the many studies given by the Russian 
school which can be found for instance in [9]. Some non-homogeneous examples are 
found in m- 

For other Hamiltonian systems on non semisimple Lie groups only few examples 
and generalizations are known in the case of the geodesic flow. This is the situation 
for nilpotent and solvable Lie groups or even their compact quotients which are 
locally homogeneous manifolds. 
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On the one hand Butler proved the integrability of the geodesic flow on Heisenberg- 
Reiter 2-step nilpotent Lie groups, among them one can hnd the Heisenberg Lie 
group [4]. He also found an algebraic condition for non-integrability. On the other 
hand, Eberlein started a study of the geometry concerning the geodesic flow on 
2-step Lie groups following his own and longer study in this topic, giving a good 
material and references in Ell]. This study of Eberlein is much more general and 
is mixed with many other geometrical questions. 

In the present paper we concentrate in the geodesic flow of Lie groups endowed 
with a metric invariant by left-translations. In the hrst part we write the basic 
definitions and get general conditions and formulas for the involution of first integrals 
making use of the Lie theory tools, that is, assuming some natural identifications 
which are not present in the work of Butler. However, all this stuff is good explained 
along the work. We write the underlying results on the tangent Lie groups and put 
special emphasis on 2-step nilpotent Lie groups, we take as nilmanifold after [22], 
for which there exists a developed geometrical theory and several examples and 
applications, see [6|. 

We apply the results we get to the case of the Heisenberg Lie group H„ of di¬ 
mension 2n -|- 1 equipped with the canonical left-invariant metric. Although this is 
a naturally reductive space the methods of Thimm do not apply in this case. 

One of our motivations is to investigate the nature of the first integrals one can 
construct. We prove that all the hrst integrals we get can be visualized on the 
isometry group. In fact, this is the case of quadratic polynomials which are invariant, 
so as hrst integrals arising from Killing vector helds. Recall that, given a Killing 
vector held X* on a Riemannian manifold M, one has a hrst integral on the tangent 
Lie bundle TM dehned by fx*{v) = {X*,v). We proved that 


(i) There is a bijection between the set of quadratic hrst integrals of the geodesic 
how on H„ —with the canonical metric— and the Lie subalgebra of skew- 
symmetric derivations of the Heisenberg Lie algebra so that involution of 
quadratic hrst integrals would correspond to a torus of skew-symmetric deriva¬ 
tions (Theorem 4.1). Actually, a general formulation of quadratic polynomials 
on a 2-step nilpotent Lie algebras to be hrst integrals is found so as the pairwise 
commutativity condition. 

(ii) The linear morphism X* fx* builds a Lie algebra isomorphism onto its 
image. 


Item (ii) gives an answer to the question formulated in |2T] and it is the hrst 
example we found of this situation. 

Making use of all results we found a new family of hrst integrals for the geodesic 
how on H„. 

After that, we consider an arbitrary lattice A and, passing to the quotient, we 
explain the integrability of the geodesic how on the compact space A\H„. Einally, 
we consider any left-invariant metric on the Heisenberg Lie group H„ and prove that 
the corresponding geodesic how is also completely integrable. 
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2. The canonical symplectic structure 

Let M denote a differentiable manifold. Its cotangent bundle T*M admits a 
canonical symplectic form 12 constructed as follows: for every ^ G Tr^{T*M) define 
the 1-form 0 by 

0(C) := i?(d7f(C)), 

where If : r*M —)■ M is the canonical projection which assigns to rj £T*M the base 
point p G M. Now the symplectic form of T*M is given by 

Q := —d&. 

If M is a Riemannian manifold with metric tensor (•,•), one induces a symplectic 
form Q in the tangent bundle TM of M. In fact, take 

0(C) := (u,d7r(C)), 

where C £ T^iTM) and tt: TM ^ M denotes the canonical projection, and then 
define 

n := —dQ. 

It is easy to see that both 0 and are the pull-back by the natural diffeomorphism 
(/) ; TM —)■ T*M given by (/>(u) = (u, •), that is 0 = </>*(©) and Q = (/)*(fi). 

We shall concentrate the study in the case where M is a Lie group endowed with a 
left-invariant metric. Let N denote a Lie group endowed with a left-invariant metric 
(•,•). Let n be the Lie algebra of N. The tangent bundle TN ot N is parallelizable 
and it is identified with 

TAT ~ AT X n = {(p, T) ; p G AT, y G n}. 

With this identification, we consider the pair (p, Y) as the tangent vector dLp{Ye) G 
TpN. Recall that TA^~A^xnisa Lie group regarded as the direct product of N 
and the abelian group n and so 

Tip,Y)iTN) ^ n X n = mV) : U,V e n}. 

When the base point is not clear from the context, we will use the notation 
{U,V)(jy^Y) foi' a tangent vector in T(py)(TN'), p ^ N,Y G n. The pair {U,V) 
denotes a left-invariant field on TN (i.e., the left-invariant vector field U on N in 
the first coordinate, and the constant field R in n in the second one). 

With this identification, and since the metric on N is left-invariant, the canonical 
1-form on TN has the following form 

eip,Y){u,v) = {Y,u). 

Recall that the curve on TN starting at (p,Y) 

c{t) = (pexp(t[/),ty -|- y) 
has initial speed (t/, V) G T(j,y){TN). So, 

{u,v\p^Y)&{u',v') = ^ = I ^{tv + Y,u'), 

and thus 

{u,v)(^p^Y)e{u',v') = {v,u'). 
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Finally, 

V), ([/', V')) = -dQiiu, y), ([/', V')) 

= -{([/, y)(p,y)©([/',F0 - {U',r){p,YMU,V) 

= -(F, [/') + {U, V') + 0(p,y)([^, U'],0), 

and therefore 

(2.1) fl(p,y)(([/,F),([/',F')) = {U,V') - {V,U') + {Y, [U,U']). 

Remark 2.1. A similar work can be done in the cotangent bundle 

T*N ~ iV X n* = {(h, a) : /i G IV, a G n*}. 

Here the pair {h, a) is identified with aodL^-i l/i, after identifying n ~ T^N as usual. 
Notice that if we move h but fix a we get a left-invariant 1-form on N. Using similar 
notations and identifications as in the case of the tangent bundle, one can get: 

Tih,a)iT*N) = n X n* = {(U,/3) : U G n, /3 G n*}, 

and hence, 

eih,a){U,^)=CXei6Lf,-lUh). 

A curve in T*N with initial speed {U,/3) G T(/j q)(T*A 1) is given by 

c{t) = {hexp{tU),tf3 + a). 

Analogously to the tangent bundle one can show that 

and 

n(^h,a)m(3),{U',(3')) = -/?([/') +/3'(U) + «([[/,[/']). 

The Lie group N is equipped with the product metric of x n which is invariant 
under the left-translations Lg for g £ N. In other words the Lie group N acts on 
TN hy isometries. From now we shall assume that TN is endowed with this metric. 
The next definitions introduce the Hamiltonian vector field and the gradient field. 

Definition 2.2. Let / : TiV —)■ R be a smooth function. 

• The Hamiltonian vector field Aj of / is defined by 

(2.2) d/(p,v)(W^) = ^{p,Y)i^f, for W G TpN, 

with Hamiltonian equation 

(2.3) c'{t) = Xf{c{t)) for a smooth curve c : RTA. 

• The gradient field for /, denoted grad/, is the vector field on TA given by 

(2.4) d/(p,y) (IT) = (grad(p,y) /, IT) for IT G T^A. 
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Let grad(p y) / = {U,V), for U,V G n denote the gradient field of a smooth 
function / : TAi —)• R. We shall compute the Hamiltonian vector field in terms 
of the gradient vector field. 

Let us denote Xf{p, Y) = {U", V”). If (17', V') is an arbitrary element of T(^py)TN 
we have that 


(H, U') + (17, V) = (grad(p,y) /, ([/', V')) 

= df^py){U',V') 

= n^p^Y){{u",v"),iu',v')) 

= {U",V')-{V",U') + {Y,[U",U']). 

Now choose U' = 0, so that the equalities above reduce to 

{V,V') = {U",V') 

implying that U” = V. By making use of this in the relations above one gets 

(2.5) {U, U') = -{V'\ U') + (y, [y, U']) for all U' G n. 

Hence the Hamiltonian vector field for / : TA^ —)• R is given by 

(2.6) Xf{p,Y) = {V,-U + ^d\V)Y) for grad(p,^) / = (H, H) 

where ad*(y) denotes the transpose of ad(y) relative to the metric on n. 

Recall that a symplectic structure on a symplectic manifold (A^, H) defines a Pois¬ 
son bracket on C°°{M) in the following way. Let /, s' : TAi —)■ R be smooth functions, 
then its Poisson bracket is given by 

(2.7) {/,s}(p,y)=%,y)(X^,X5) 

and one says that / and g are in involution if {/, g} = 0, alternatively /, g Poisson 
commute. 

Lemma 2.3. Let {N, (, )) denote a Lie group equipped with a left-invariant metric 
and let / : TAf —)■ R denote a smooth function. If one denotes the gradient of f by 
grad(p^y) f = {U,V), then 

the Hamiltonian vector field of f is given by 

Xf{p,Y) = {V,ad\V){Y)-U) 

where ad*(y) denotes the transpose of ad(y) with respect to the given left- 
invariant metric on N. 

The Poisson bracket of a pair of functions f,g:TN^Wis given by 

{f,9}iP,y) = (A grad/, A grad S') 

= {V',U)-{V,U>)-{Y,[V,V']), 


(^) 

( 2 . 8 ) 

(ii) 

(2.9) 


where A{U,V) = {V,U). 
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Proof. The statement (i) was proved above. For the proof of the second statement 
assume grad(-p y) / = {U, V) and grad^p y) g = {U', V). Then 


{f,9]{p,Y) 


f ^ Yg) 

0(p,y)(F,ad‘(F)(y) - [/), {V',^d\V'){Y) - U') 
{V',U)-{V,U') + {Y,[V’,V]) 
-{{V,U’)-{V',U) + {Y,[V,V'])) 


which concludes the proof. 


□ 


2.1. Example: 2-step nilpotent Lie groups. Assume is a connected 2-step 
nilpotent Lie group with Lie algebra n. Whenever N is furnished with a left-invariant 
metric (, ), one can read several geometrical features of N at the Lie algebraic level. 
In fact, the Riemannian metric is determined at the identity element so that the Lie 
algebra n can be decomposed into a orthogonal direct sum 

n = 0 0 3 with 0 = j-*- 

where as usual 3 denotes the center of n. The Lie bracket on n induces for Z E 3 the 
skew-symmetric linear map j{Z) : d —)• d given by 

(2.10) {[U^ViZ) = {j{Z)U,V) for ZE 3 ,C/,FEd. 

Conversely, let (a, (, )a) and (b, (, )[,) denote vector spaces endowed with respec¬ 
tive inner products. Let n = b © 0 denote the direct sum as vector spaces and let 
(, ) denote the product metric on n given by 

(2-11) (.)u. = (.>. (.>L.. = (.>b {o.l')=0. 

Let j : 0 —?■ End(b) be a linear map such that j{Z) is skew-symmetric with respect 
to (, )b for every Z E a. Then n becomes a 2-step nilpotent Lie algebra if one defines 
a Lie bracket by the relation in (I2.10h and so that a is contained in the center of n, 
0 C 3 ; actually 

3 = 0 0 {F E b : [R, [/] = 0 for all C/ E b} 

= a© Czenker j(Z). 

By using the translations on the left, the corresponding connected Lie group N is 
endowed with a Riemannian metric. 

The inner product (, ) produces a decomposition of the center of the Lie algebra 
n as a orthogonal direct sum as vector spaces 

3 = kerj ©(©(n), 

where C'(n) denotes the commutator of n and the linear map j is injective if and 
only if there is no Euclidean factor in the De Rahm decomposition of the simply 
connected Lie group {N, (, )) (see [iQj). 

Let n be a 2-step nilpotent Lie algebra. We say that n is non-singular if ad(A) : 
n ^ 3 is surjective for all X E n — 3 . One can show that the following properties are 
equivalent - see for instance [6]: 

(i) n is non-singular; 

(ii) for every inner product (•, •) on n and every nonzero element Z of 3 the linear 
map j{Z) is non-singular; 
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(iii) for some inner product (, ) on n and every non-zero element Z of ^ the linear 
map j(Z) is non-singular. 

Let us denote by VLt, (resp. W^) the o-component (resp. 3 -component) of an ele¬ 
ment VL E n. Hence for arbitrary elements Y, V,U' £ n one has 

{Y,[V,W]) = {Y„[V,,W,]) = {j{Y,)V„W,), 

and the symplectic structure on TA^ can be written in the following way 

n^p,Y)mv),iu',v')) = {u,v') - {v,u') + {jiY,)u„u'). 

In particular if N is 2-step nilpotent and / : TA^ —>■ M is a smooth map with 
gradient field grad(p y) / = {U, V) then Lemma [32] says 

(2.12) Xf{p,Y) = {V,jiY,)V,-U). 

Furthermore, for f,g : TN —>■ M smooth maps with respective gradients given as 
grad(p y) / = {U, V) and grad(p y) g = {U', V'), the Poisson bracket is 

(2.13) {f,g]{p,Y) = {V',U) - {V,U') + {j{Y,)V;,V,). 

3. The geodesic flow 

Let N denote a Lie group with tangent bundle TN. In the following paragraphs 
we will denote by (•, •) both the left-invariant metric on N and the product metric 
on TA^ ~ A^ X n, that is, the product of the left-invariant metric on N and the inner 
product on n induced by the left-invariant metric of N. 

The energy function E : TA^ —)• M is defined by 

E{p,Y) :=^{Y,Y)p = ^{Y,Y), 

since (•, •) is a left-invariant metric. 

Notice that 

dE(p,y)([/,E) = fX'^{peMtU),tV + Y) 

(3-1) = m{tHV,V) + 2t{Y,V) + {Y,Y)) 

= {y,v). 

The gradient field of the energy function grad E is the vector field on TAl implicitly 
defined by the equation 

dE(p,y)(t/, E) = (grad(p,y) E, {U, V)), 

and thus the computations given in ()3.1I) imply 

grad(p^y)E = (0 ,y). 

The geodesic field on TA^ is the Hamiltonian vector field of the energy function 
Xy- The Equality ()2.8p for the energy function implies 

XE(p,y) = (T,ad*(y)y), 
and whenever N is 2-step nilpotent we get 
(3.2) XY{p,Y) = {Y,j{Y,)Y,). 

Let (M, (, )) denote a complete Riemannian manifold with tangent bundle TM. 
Eor each v £ TM and t £ M. define r*(u) = 7 (,(t) the velocity at time t of the 
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unique geodesic with initial velocity v. The fact that M is complete implies that the 
geodesics of M are defined on M. Thus for every v G TpM the curve on TM given by 
r*(u) is defined for all t G M. One can check that T* o T'^ = for all s, f G M. The 
geodesic vector field is taken as the vector field on TM with flow transformations 

{T^. 

On the other hand one has the flow of the Hamiltonian vector field Xe in TM 
which is determined by the energy function E : TM — )• M, that is E{v) = ^ < u, u > 
for all V G TM. Now the flow of Xe coincides with {T*}. See for instance Section 5 

in [8]. 

This also applies on a Lie group N endowed with a left-invariant Riemannian 
metric. Moreover in EE] one can see the definition and properties of the geodesic 
flow in the Lie algebra n and the relationship between this with the geodesic flow 
defined above in terms of the Gauss map G : TN ^ n. 

Let (M, (, )) be a Riemannian manifold. The so called first integrals of the geo¬ 
desic flow are the functions f : TM ^ which Poisson commute with the energy 
function. This gives 

0 = {/, E}{v) = dUXE) = XeU) = j-J o r (u), 

from which it is clear that first integrals are functions which are constant along 
geodesics. 

First integrals for the geodesic flow can be constructed from Killing vector fields. 
In fact for a Killing vector field X* on M, the function fx* : TM ^ M given by 


(3.3) 


fx*{v) = {X*,v) 


becomes a first integral of the geodesic flow. In fact if 7 (f) is a geodesic on M and 
V denotes the Levi Civita connection, then 


L(x-(7(t)),y{t)) = {Vy(„jf,7'(t)) + /x-(7(t)),Ly(t)) 

= (Vy,„X-,7'(()) =0 


since is skew-symmetric. So {{X*{'y{t)),^'{t)) is constant, and this implies 

that fx* (v) = {X{7r{v)),v), is constant along the integral curves of the geodesic flow. 


Lemma 3.1. Let M be a Riemannian manifold and let X* be a Killing field on M. 
Then fx* is a first integral of the geodesie flow. 


In the setting of the previous section, let N denote a Lie group and let / : TX —?• M 
be a smooth function with gradient field grad(^pY){f) = {U,V). In view of the 
formulas in Lemma 12.31 the function f is a first integral of the geodesic flow if and 
only if 

H(p,y)(XE,X/) = n^py){{Y,ad\Y)Y), {V, -U + ad*(K)y)) 

(K ,-U + ad\V)Y) - {ad\Y)Y, V) + (K, [K, K]) 

(y, -U + ad\V)Y) - (ad*(y)y, V) + (ad‘(y)y, V) 

{Y,-U + ad\V)Y). 

This proves the following result. 


0 = 

(3.4) = 
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Lemma 3.2. Let {N, {, )) denote a Lie group equipped with a left-invariant metric. 
Then the smooth function f : TN with gradient 

grad(p_y)/ = {U,V) 

is a first integral of the geodesic flow if and only if for all {p, Y) &TN it holds 
(3.5) {Y,U) = {Y,[V,Y]). 

In particular if N is 2-step nilpotent, then f is a first integral if and only if 

{Y,U) = {j{Y,)V,,Y,). 

We shall say that / : TA^ —)• M is invariant if /(p, Y) = f{g-p, Y) for all g,p G N, 
Y G n, that is, / is invariant under the left-action of N into TN. 

We make use of the Equation ()3.5p to find invariant first integrals on Lie groups. 

Proposition 3.3. Let {N, (, )) be a Lie group endowed with a left-invariant metric 
and let TN denote its tangent bundle with the product metric. 

• Let fzo : TN be defined by 

fzo{p,Y) = {Y,Zo). 

Then fzo is a first integral of the geodesic flow for all Zq G 3 . Moreover, the 
family {fzi '. Zi G 1 } is a commutative family of first integrals. 

• Let A : n ^ n be a symmetric endomorphism of n and let 

gA{p,Y) = ^{Y,AY). 

Then qa is a first integral of the geodesic flow if and only if 

Q = {Y,[AY,Y]). 

Proof. It is easy to see that grad^p = (0,Zo) and so the Hamiltonian vector 

field is = (Zo,0). It follows from Equality ()3.5p that fzg is a first integral of 
the geodesic flow. Moreover, 

= ^(p,y)((^o,0), (Zi,0)) = (y, [Zo,Zi]) =0, 

which completes the proof of the hrst item. 

We will deal now with quadratic first integrals of the geodesic flow. 

Let A : n —)• n denote a symmetric endomorphism of n and let dehne 

gAip,Y) = ^{Y,AY). 

An elementary calculation gives 

dgA\ip,Y){U,V) = {AY,V) 

and hence 

gvad^ y){9a) = (0,Ay). 

Then by Equation (|3.5p . we have that qa is a first integral of the geodesic flow if 
and only if 


o = (y,[Ay,y]). □ 
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Remark 3.4. The last equality says that for A symmetric on a 2-step nilpotent 
Lie algebra, the Poisson commutativity of gA with the energy function depends 
essentially on the restriction of ^ to 0 . 

Theorem 3.5. Let {N,{, )) denote a 2-step nilpotent Lie group endowed with a 
left-invariant metric and let TN denote its tangent bundle with the product metric. 
Let n be the Lie algebra of N. Let ^ : d —)• d denote a symmetric map with respect 
to {, )o and let gA ■ TN ^ R be given by 

gA{p,Y) = ^{Y,AY), 

where AZ = 0 for all Z ^ i. Then 

(i) gA is a first integral of the geodesic flow if and only if for any basis {Zi ,..., Zm} 
of I one has 

[J(Zj ),^]=0 for all i = 1,... ,m. 

(a) Let gA and gs be a pair of first integrals, then they Poisson commute if and 
only if 

j{Zi)AB = j[Zi)BA for all i = 1,... m. 

In particular if n is non-singular, then 

{9A,gB} = 0 if and only if [A, B] = 0. 

Proof, (i) Proposition 13.31 says that gA is a first integral if and only if for T G n one 
has 0 = {j(YflAYx,,Y„) which for a basis {Z^} of 3 is equivalent to 

(3.6) 0 = {j{Z,)AY„,Y„). 

In fact for the proof of (^) take Y = V + Zi with P G d to get (|3.6p . The converse 
is clear since j{Yfl = aij{Zi) for a* G M and {Zi, ..., Zra\ a basis of 3 . 

Now take = P -|- P' G d and write 

0 = (j(Z,)A(P + P'), P + P') = (j(^*)4lP, P') - {Af{Zf)V, V) 

which implies that j{Zi)A = Aj{Zi) for all i = 1,... m. 

(ii) Assume that gA-,9B is a pair of first integrals of the geodesic flow (so that 
they satisfy (i) above). Then Lemma 12.31 savs that 

0 = {gA, 9 B] = {Y,[AY,BY]) 

= {j{YflAY,,BY,) 

= {Bj{Y,)AY„Y„) 

Firstly as in the proof of (i) the last equality is equivalent to 0 = {Bj{Zi)AY^,Yf) 
for alH = 1,... m. Also as above Bj{Zi)A must be skew symmetric so that 

{Bj{Zi)Af = -Aj{Zi)B = -Bj{Zi)A for alH = 1,... ,m. 

Since the functions gA,9B are first integrals one has 

j{Zi)AB = j{Zi)BA 

which finally implies the statement. The last assertion follows from the last equality. 

□ 
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Let be a Lie group of dimension n. We shall say that the geodesic flow is 
completely integrable if there exists n smooth functions fi, f 2 , ■ ■ ■, fn ■ TN —)• M 
such that 

• the gradients of fi, f 2 , ■ ■ ■, fn are linearly independent on an open dense 
subset of TA^; 

• 0 = {fi, fk] = {fi,E} for alli,k = l,...,n. 

We now consider a sufficient non-integrability criterion by Butler [5]: 

Definition 3.6. Let n be a 2-step nilpotent Lie algebra. 

(i) For A € n*, let ua := {X G n/ad*(X)A = 0}. 

(ii) A A G n* is called regular if ua has minimal dimension. 

(iii) n is called non-integrable if there exists a dense open subset W of n* x n* 
such that for each (A,/i) G W, both A and fi are regular and [nA,n^] has 
positive dimension. 

In view of the metric on n any element A G n* can be realized in the form 
X = {V + Z,-) for unique IL G 0 and .Z' G 3 . By identifying A G n* ^—5- V Z ^ n 
we get 

nv+z = {A: Gn : (F + Z,ad(X)-) =0} 

= {X en : {V + Z, ad{X){V' + Z')) = 0 for all V + Z' e n} 

= {X en :{Z, [X, V']) = 0 for all V G o}, 

from which it is clear that 3 C tta for all A G n*. Moreover if we denote X = V + Z 
then 

B + Z G nv+z ^ {j{Z)V, •) = 0 4 ^ {V,j{Z)V') = 0 for all V' G 0 . 

Hence if n is non-singular then V = 0 since Imj(Z) = t). In this situation ny+z 
has minimal dimension for Z 7 ^ 0 and [nA,n^] = 0. One gets the following result. 

Lemma 3.7. Let n be a non-singular 2-step nilpotent Lie algebra, then n cannot be 
a non-integrable Lie algebra. 

Thus non-singular Lie algebras are candidate to have a completely integrable 
geodesic flow in view of the next result. 

Theorem ([5], Theorem 1.3). Let n be a non-integrable 2-step nilpotent Lie algebra 
with associated simply connected Lie group N. Assume that there exists a discrete 
cocompact subgroup A of N. Then for any such A and any left-invariant metric g 
on N, the geodesic flow of {X\N, g) is not completely integrable. 

An example of a non-singular Lie algebra is the Heisenberg Lie algebra t)n cor¬ 
responding to the Heisenberg Lie group H„. We shall see in the next section that 
the geodesic flow on any quotient A\H„ is completeley integrable for any induced 
left-invariant metric. 

Remark 3.8. Topological obstructions to integrability on two-dimensional surfaces 
were given in |12l 114] . In any dimension on non-simply-connected manifolds by 
Taimanov m- 
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4. On THE FIRST INTEGRALS OF THE GEODESIC FLOW ON THE HEISENBERG 

MANIFOLDS 


A Riemannian Heisenberg manifold is a compact manifold given as a quotient 
A\Hri, where Hn denotes the Heisenberg Lie group of dimension 2n + 1 and A < H^ 
is a discrete cocompact subgroup. This notation coincides with that one in 

The Heisenberg Lie group H„ is the simply connected Lie group constructed over 
the smooth space together with the multiplication map: 

(u, z){v', z') = (v + v\z + z'- V, v' G ^ G M, 


where J is the real linear map representing the multiplication by i once we identify 
with C"". It is clear that a basis of left-invariant vector fields is given for 
A: = 1,... ,n by 

Zi{p)=dz, A2fc_l(p) = 5x2fc_i - X2k{p) = dx2t^ + ]^P2k-ldz-, 

where Xi are the canonical coordinates on and with dy we denote the partial 
derivative So any left-invariant vector field U on H„ has the form 

2n 2n / ^ n \ 

U{p) = aiXi + cZi = '^ Qidxi + ( c + 2 '^{-o-2i-iP2i + a2iP2i-i) | d^. 

i=i i=i V i=i J 

Note that Yll=i{~^ 2 i-iP 2 i + CL 2 iP 2 i-i) = A) where uj is the canonical sym- 
plectic form of and P = {pi,p2 ,... ,P 2 n), A = (oi, 02 ,..., a 2 n)- 

Let (, ) be the inner product on ()„ which turns the set {Xi, X 2 , ■ ■ ■, X 2 m Zi} 
into an orthonormal basis. In canonical coordinates of we get the following 

metric: 


2n ^ ^ 2n 

dz^ + dzdxk + (1 + jvDdxl] + 2 X] i-^)^^''ykysdxkdxs 

k=l k^s,l 


for t = 21 — 1 
for t = 21. 

We also denote by (•, •) the left-invariant metric on H„. Then the restriction of 
(•, •) to the Lie algebra \]n gives a decomposition 


where yt = 


X 21 

X 21-1 


fin = 0 ® 3 as orthogonal sum of vector subspaces, 


where 0 is the subspace spanned by Wi, A 2 ,... ,X 2 n and 3 = MZi is the center of 
f}n. For a generic left-invariant vector field U denote by U„ and the respective 
projections onto d and 3 . One can see that the integral curve of U at the point 
P= {Pv,Pz) is given by 

1p{s) = {Pv + sU„,Pz + sU^ - ^splJUr,) 

via natural identifications at the corresponding vector spaces. 

It is the clear that the exponential map exp : f)„ —)• H„ is given by 

exp([/) = U 
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via identifications: at the left-side U is given in terms of a left-invariant basis, while 
at the right side the vector U E is given in canonical coordinates. 

Under the Lie bracket on X(H„) the set of left-invariant vector fields give rise to the 
Lie algebra called the Heisenberg Lie algebra whose basis {Xi,X 2 , ■ ■ ■, X 2 n, Zi} 
satisfies the non-zero Lie bracket relations 

[X2i-i,X2i\ = Zi, i = 

so that for Z ^ i, the skew-symmetric endomorphism j{Z) E so(t)) as in (12.1011 takes 
the form 

j(Z) = (Z,Zi)J 

where J := j{Zi) has the following matrix 

/O -1 \ 

1 0 

0 -1 

V 10/ 

in the basis {Xi,X 2 , ■ ■ ■, X 2 n} of 0. 

Since the dimension of the center of i)n satisfies dim 3 (f}„) = 1, we can only define 
one independent linear first integral for the geodesic flow: 

fzAp,y) = {y,zi). 

As we said before, we have that 

grad(p_y) fzi = (0,Zi), Xf^^{p,Y) = (Zi,0). 

As known, the isometry group of a simply connected nilmanifold M is the semidi- 
rect product Iso{M) = K t< N where N is the nilradical of the isometry group which 
acts simply and transitive on M. Naturally M can be identified with N so that the 
homogeneous metric on M turns into a left-invariant metric on N. Moreover the 
isotropy subgroup K consists of isometric automorphisms. See [22] for more details. 

Let K be the subgroup of isometric automorphisms of the Heisenberg Lie group 
H„, its Lie algebra i is given by 

(4.1) t = {H E 5o(d, (, )d) such that [J, H]=0}, 

see for instance m- 

It is clear that the action of H„ on itself by isometries is given by the translations 
on the left. 

Let us search for the quadratic homogeneous polynomials which are first integrals 
of the geodesic flow. For a symmetric map A : d —)■ d the quadratic polynomial 
given by 

gA{p,y) = {AY,Y) 

is a first integral of the geodesic flow on the Heisenberg Lie group H„ if 

[J, A] = 0 by Theorem 13.51 
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Take ^ : D —)• D a symmetric map such that [ J, A] = 0 then the map given as 
B := JA E so(t)(, )d) also satishes 

[J,B]=Q 

so that B belongs to the isometry algebra of In fact, as explained above 

/so(H„) = K tK H„ where K denotes the isotropy subgroup consisting of isomet¬ 
ric automorphisms. Its Lie algebra is iso(H„) = fi © [)n where t is the Lie algebra of 
K and it consists of skew-symmetric derivations, see (SI]). 

Conversely take B G t and define a symmetric map ^ : d —tJ by ^ := JB. Then 
gA is a first integral of the geodesic flow. 

Theorem 4.1. There is a bijeetion between the set of quadratic first integrals of 
the geodesic flow and the Lie subalgebra of skew-symmetric derivations of given 
by 

(4.2) A —^ fi(A) := JA 

so that {fi'Aiifi'Aa} =0 if and only if V'(^ 2 )] = 0 in so(d, (, )„). 

Remark 4.2. Note that is not a Lie algebra morphism. 

Remark 4.3. The Gauss map G : TH„ —?■ is an anti Poisson map. The metric 
(, ) and the Lie bracket on f)„ define a Poisson structure on the Heisenberg Lie 
algebra t)n- For the energy function on f}„ given by E{v) = \{v,v) one also has 
its Hamiltonian vector field called the geodesic vector field, which corresponds in a 
natural way to the geodesic flow means of the Gauss map G ; TH„ —)• f)n- In fact 
if / is an invariant function on TH„ which corresponds to a function / : f)„ —>■ M 
then G{Xf{p,v)) = —Xj{v) where Xj denotes the Hamiltonian vector field for / 
while Xj for /. If P* denotes the geodesic flow for Xe and T* the flow of Xg then 
G o T-* = f * o G for all t E M. See for more details. A formula for geodesics on 
the Heisenberg Lie group can be found in [6]. 

The coadjoint orbits on the Heisenberg Lie algebra were described in |16] were 
making use of the Adler-Kostant-Symes scheme they were used for a study of the 
motion of n harmonic oscillators near an equilibrium position. 

From the result above it is clear that to find the symmetric maps S' : d —>• d such 
that [J, S] = 0 is equivalent to determine the skew-symmetric maps T E so(d) such 
that [J,T] = 0. Choose Ai, A 3 ,..., A 2 n-i, A 2 , A 4 ,..., A 2 n the ordered basis of d 
and write T and J in this basis. A general matrix T E 5 o( 2 n) of the form 

f Qi Q 2 \ 

\-Ql Q3) 

with Qi,Q 3 E 50(n) real submatrices of T, commute with J 



[J, T] = 0 if and only if Qi E so(n), Qi = Q 3 and Q 2 = Q 2 - 
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Take the quadratic first integrals given by the symmetric linear maps : 0 —d 
which correspond to 

(4.3) AiY = {Y,X2i-i)X2i-i + {Y,X2i)X2i, i>l. 

It is not hard to prove that [J,= 0 and [Ai^Aj] = 0 for all z,j = 1,... ,n, 
which gives { 51 ^. T) = 0 by Theorem 13.51 

Note that the first integrals /z^, gAi (f > 0), are algebraic, that is invariant under 
the action of H„ on TH„ ~ H„ x t)n- 

Recall that Killing vector fields built the Lie algebra of the isometry group of a 
given homogeneous nilmanifold. So for the Heisenberg Lie group we have 


Iso{Rn) = K X 


where H„ is the normal subgroup of /so(H„) corresponding to the translations on 
the left by elements of H„, K is the subgroup of orthogonal automorphisms and the 
action of K on H„ is given by the evaluation map. 

Let T E so(d) such that [J,T] = 0 and pT be the automorphism of H„ given by 

Pt{p) = {e^Pv.z) 

where py is the projection of p onto and is the usual exponential map: 

• Since pT is a group homomorphism the fact that it is an isometry 
at the identity element implies that pT is an isometry of H„. 

For this isometry the corresponding Killing vector held computed as 


is given by 
(4.4) 


xKp) 


A. 

ds 


Pt{p) 


d 

ds 


{e^'^py,z) 

0 


X^{p) = TW, - ^{TW^, JWy)Zi 


where T E so(d), so that exp W = p for the exponential map exp : t)n —>■ H„ and VLo 
denotes the projection of W onto 0 with respect to the orthogonal decomposition 
Ifn = 0 ©T 

The Killing vector helds which generate as subalgebra of the isometry algebra 
corresponds to the right-invariant vector helds given by 

Zlip) = d^, X^i-i{p) = d^^._^ + ^P2id^, X^i{p) = d^^, + ^p2i-id,. 

Explicitely, for X^ with, k = 0,..., 2n, corresponding to the left-translations 
isometries, we get the next hrst integrals - see (I3.3jl : 


F 2 ^-l{p,Y) = {X 2 i-i + {W„X 2 i)Z^,Y) for X 2 V 1 , f = 1 , • • • ,n 

F 2 ^ip,Y) = {X 2 i - {W„X 2 i-i)Z^,Y) for ^ 2 *^, i = l,...,n 

which for k = 1,... ,2n can be rewritten as 

Fk{p,y) = {Y-JiY,)Wy,Xk) 

^ ’ = {Y,Xk)-{J{Y,)W„Xk), 

while for Zf one has the function 

fz,{p,y) = {Zi,y) 
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which was introduced in the previous section, see Proposition 13.31 

For the Killing vector field corresponding to an element T of the isotropy 
subgroup K C /so(H„) one has the first integral 

(4.6) Friv^y) = {TW„Y,) - VF„)(Zi,K) 

where A is the symmetric map of the bijection (I4.2jl for T = JA. 

For these functions one gets the following gradient vector fields 

grad(p,y) Ft = {-{Zi,Y)AW^ - TY„, TW^ - ^{AW,, W,)Z,) 

(4.7) grad(py)Ffc = {{Y, Zi)JXk, Xk + {W, JXk)Zi) foi k = I,... ,2n 

grad(p_y)/zi = ( 0 :- 2 ’i) 

which shows that the gradients of fzi, F^ for k = 1 ,... 2 n are linearly independent. 

The formula in (j2.13p says that for j, k G {0,..., 2n} the Poisson bracket is given 
by 

(4.8) {Fj,Fk}{p, Y) = {Y, Zi)(Xfc, JXj) 
so that the only non-trivial brackets are 

(4.9) {F2i-i,F2i} = fz,. 

Let PTa^PTi denote a pair of isometries as above with Ta = JA, Tb = JB. Recall 
that Ta,Tb G 5o{v) and 0 = [J,Ti] = [J,Tj]. Take the corresponding functions 
Ft^ , Ft^ as in (j4.6p with corresponding gradient vector fields grad(p ^y^FT, = iU,V) 
and grad(p y) Ft^ = {IF, V). Then by Lemma [23] one has 

{Ft^,Ft,}{p, y) = 0 if and only if {V\ U) - (K, U') + (T, [K', K] = 0). 

In particular for (p, Z \) one gets 

{Ft., Ft.Kp, y) = -{TbW„AW,) + (TaW^BW,) + (JTbW^TaW,) 

= -{JABW„,W„) + {JBAW^, LFt,) + JAW„) 

= -{JABW„W,) + {JBAW^,W^) + {JABW„W,), 

so that as in the proof of Theorem 13.51 {F^^. F'r. }(p. Y) = 0 implies [A,B] = 0 for 
Ta = J A,Tb = JB. 

Conversely it is not hard to see that (Fy., Ft^}(p, y) = 0 whenever [A, B] = 0 for 

all (p,y). 

Lemma 4.4. Let Ta = J A and Tb = JB denote skew symmetric derivations. Let 
Fxa , : TBn —>• M 6e the functions defined above. It holds 

{FTa,FTj{p, y) = 0 if and only if [^, B] = 0. 

Moreover, if [A,B\ = 0 then (Fy., 5 s}(p, y) = 0 whenever A,B are symmetric 
maps. 

Making use of ()4.7p the other Poisson brackets at every point {p, Y) for p = exp W, 
satisfy 

• {Fy,, 5 s}(p,y) = {JABY„Y,) for T, = JA. 
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• Let Ta = JA, Th = JB for A, B symmetric. 

{FT^,Fn}{p, Y) = {[Ta, n]W„, n) + {Y, Z,){JBAW„, W,) 

= {[Ta,n]W„,Y^) + {Y,Zi)\{{JBA-JAB)W^,W^) 

= ^Ta,ny 

Notice that — JBA) is the symmetric part of JAB. 

• {Ft, Fk}{p, Y) = Fkip, TY, + Y,), for A: = 1,, 2 n. 

• {/zi, FJ = 0 for all i = 1 ,..., n. 

Let Ai : X) ^ X) the symmetric maps as in (14.3p and let Ti = JAi the corresponding 
skew-symmetric maps. Then from the equations of the gradients above one can see 
that FJ'^,gA■ is a set of 2 n linearly independent functions on an open dense subset 
of TRn- 

Note that the set {JAi} correspond to a Cartan subalgebra of so(]R^") which in 
particular gives an abelian subalgebra of skew-symmetric derivations. 

We are now in conditions to prove the following result. 


Theorem 4.5. The geodesic flow on TH^ is completely integrable in the sense of 
Liouville. In fact the set 


g = {F}U{5Ajr=iU{FTjLi 

gives a family of commuting first integrals of the geodesic flow. 
The sets 


( 1 ) 

(2) and 


F={/zJU{ffAjr=lU{F2fc_i}Ll 


= {/zj U {<7Ajr=i U 

give two independent eommuting families of first integrals of the geodesic flow. 


Proof. Lemma 14.41 savs that 

• {FTi, Ftj } = 0 and {Ft^ , gAj} = 0 for Ti = JAi with i = 1,... ,n and 

• {9Ai,gAj} = 0 since [Ai,Aj] = 0, 

and since all of them are first integrals the fact that this is a completely integrable 
set follows from the condition of linearly independence of the set on a open dense 
subset of TRn- 

For the other families by making use of the information already given we have 
that 

• {F 2 k-i,F 2 j-i} = 0 = {F 2 k,F 2 j} for all k, j = I,... ,n and k j. 

We should compute the other Poisson brackets. 

ifz^gAMP^Y) = {Z^,j{YflAiY) + {Y,[Z^,AiY]) = Q, 

{gA,,F2j-i}{p, Y) = -{j{Y,)A,Y, X2j-i) + {Y, X2j-i]) = 0, 

{gA,,F2k}{p,Y) = -{j{YflA,Y,X2k) + {Y,[AiY,X2k])=0, 

which completes the proof of the theorem. □ 
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Remark 4.6. It turns out that the families given in [Ij are of the same type as the 
families in and above. While the first integrals Ft^ are new. 

Notice that in the family Q all the first integrals are obtained from the Lie algebra 
of the isotropy subgroup K C /so(Hji). 

Remark 4.7. In the simply connected case we are working on all the hrst integrals 
introduced in Theorem l4.5l are analytic. As we shall see later, this cannot be achieved 
for compact quotients of H„. 

The Poisson brackets of the functions {Ft, F^, fz^} for T G t, F^, k = 1,..., 2n 
give the proof of the next result. 

Theorem 4.8. The linear morphism between the Killing vector fields on and its 
image in C'°°(TH„) given by 

X* fx* 

as in ,71) builds a Lie algebra isomorphism between the isometry Lie algebra o/H„ 
and its image equipped with the Poisson bracket. 

Proof. In fact, a Killing vector field X* can be written as X* = X^ + Xf^ for 
T + 1/ G I © f}n- Moreover U{p) = V{p) + zZ^ and V{p) = SkXf.{p). Thus 

fx* {p, Y) = Ft{p, Y) + Y, SkFkfp, Y) + zfz, {p, Y) 

k 

SO that X* fx* is a Lie algebra homomorphism. It is injective: in fact assume 
(4.10) 0 = Ft{p, Y) + Y^ SkFkip, Y) + zfzfip, K) for all {p, Y). 

k 

By taking p = exp Zi, Y = Zi one gets z = 0. 

Take T = Zi so that Ft{p,Y) = and Fk{p,Zi) = -{JW^,Xk). 

Thus equality (14.101) becomes 

0 = -]^{AW„Wf) -Y,Sk{JW^,Xk). 

k 

If {AXj,Xj) = 0 for all j = 1,... ,2n then A = 0 and so Ft = 0. Choosing 
W = Xi one gets 0 = Y2k Sk{JXi,Xk) = ±s; which says s/ = 0 for every 1. 

If {AXj,Xj) 0 then there is j such that a^j = {AXj,Xj) 0. Take W = tXj 
for t G M. Then (|4.10p becomes 0 = —\ajjt^ — = t{—^tajj — ±Sj) which should 

holds for every t. The sign ± depends on the parity of j. Thus ajj = 0 and Sj = 0 
for all j. And X* fx* is injective onto its image. □ 

So far we know this is the first example among nilpotent Lie groups of this iso¬ 
morphism. The question appeared in [ 2 T]. It could be interesting to know if there 
are more examples of this situation among nilpotent Lie groups and if is there a 
relationship with the complete integrability of the geodesic flow. 

Remark 4.9. The metric considered in this section make of a naturally reductive 
Riemannian space. That means that there exists a Lie group of isometries G acting 
transitively on N such that its Lie algebra 0 splits as 0 = 1} © m where f) represents 
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the Lie algebra of the isotropy subgroup and m is a Ad(Lf)-invariant complement 
such that 

(x, [y, z]m) + {y, [x, z]m) = 0 for all x,y,z G m. 

The Heisenberg Lie group does not correspond to the examples in [21j . 

4.1. Riemannian Heisenberg manifolds. Let the (2?T-+l)-dimensional Heisen¬ 
berg group endowed with the standard Riemannian metric, which is defined in the 
beginning of this section. For each n-tuple r = (ri,...,r„) E (Z+)’^ such that 
I r 2 I • • • I , we define 

(4.11) = {{v,z) : V = (x,y) with x E rZ”, y E 2Z"', z E Z}, 

where x = (xi,..., x„) E rZ” means that Xj E r^Z for all z = 1,..., n. 

It follows from m that this family classifies the cocompact discrete subgroups of 
H„ up to isomorphism. This is explained in the following remark. 


Remark 4.10. In several works the Heisenberg Lie group is defined started with the 
triangular ( 2 n -|- 2 ) x ( 2 n -|- 2 ) real matrices of the form 

lx t\ 

0 In , 

0 0 1 / 

where x,y E M” and In is the n x n identity matrix. The map dehned as 
<I>( 7 (x, y, t)) = (x, y,t — ^x ■ y) gives an isomorphism with the Heisenberg Lie group 
defined above in the first paragraphs of this section. In [n] it is proved that any 
lattice on the Heisenberg group, with this matrix presentation, is isomorphic to one 
of the form 

Tr = {■y{x,y,t) : x E rZ"", y E Z", t E Z}, 

where r = (ri,..., r„) E (Z"*")” is such that ri | • • • | r„. Moreover, T^ is isomorphic 
to Ts if and only if r = s. Recall that, 4>“^(Ar) C T^. This contention is strict, 
but one can still prove that A,, is isomorphic to T^. For the sake of completeness we 
include the proof of this fact in the following remark. 


lix,y,t) = 


Remark 4.11. Recall that the center Z(Tr) of F^ is { 7 ( 0 ,0,t) : t E Z} for all r. So 
Z{Tr) is cyclic with two generators 7 ( 0 , 0 , 1 ) and 7 ( 0 , 0 ,- 1 ). As it follows from [TT] . 
for each i = 1 ,..., n we have that 

7 (rjej, 0 , 0 ) 7 ( 0 , e*, 0 ) 7 (rjei, 0 , 0 )“^ 7 ( 0 , e*, 0 )“^ = 7 ( 0 , 0 , 1 )'’*, 

and for all s E (Z'*')"' such that si | • • • | Sn, s r, there are no elements 
71,72,..., 72n E F^ such that 

72i-i72i72“iii72'i^ = 7(0,0,±1)’'L 

This shows, in particular, that F^ is not isomorphic to F^ if r 7 ^ s. Now, from the 
classification theorem in we have that A^ is isomorphic to F^ for some s. But 

(riej,0,0)(0,2ei,0)(-riei,0,0)(0,-2ei,0) = (0,0, r*) = (0,0,1)'’S 


where (0, 0,1) is one of the two generators of Z(Ar). Therefore, A,, is isomorphic to 

F,. 
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Let be as defined in ()4.1ip and let assume that H„, Ar\H„ are endowed with 
the standard metric. Since the quotient projection vr : H„ —Ar\Hn is a Riemannian 
submersion and furthermore a local isometry, we can identify the tangent bundle of 
A\H„ with (A^\Hn,) x tin- The projection vr maps geodesics into geodesics and the 
energy function E of T(Ar\H„) is related to the energy function E of rH„ by 

E{Kp,Y) = E{p,Y) = ]^{Y,Y) 

and this does not depend on the given representative. 

Since the integrals fz^ , QAi of the geodesic flow of TH^, as given in Theorem 14.51 
do not depend on the first coordinate, they descend to first integrals 

fzd^rP,Y) = fzi{p,Y), gA,{Kp,Y) = gAi{p,Y) 
of the geodesic flow of T(Aj.\H„). Moreover, such first integrals are in involution, 
since for all /, 5 G C'°°(r(Ar\H„)) we have 

{/o7r,5 0 7r} = {f,g}o7r. 

Note that the integrals F^, k = 1,..., 2n, from Theorem 14.5l do not descend to the 
quotient. However, let {p, Y) G TH„ and q G A^, say p = (x, y, z) and q = (x', y', z'), 
with x,y £ M"", z G M, x' G rZ"", y' G 2Z"' and z' G Z. Take W,W',W'' G f}n such 
that exp W = p, exp W' = q and exp W" = qp. Observe that W” = Wx, + W^. Thus 
one has 

Fk{qp,Y) = (y,Xfc) - (y,Zi)(J(TT„' + fT„),Xfc) 

= Fk{p,Y)-fzAP,Y){JW',Xk). 

Since {.JW^,Xk) £ h we have that 

Fk {qp, Y) = Fk {p, Y) mod fz, {p, T)Z 
and since fz, is a first integral of the geodesic flow we have that the function 

'■■-"“-(-IS) 

descends to Ar\H„ and is constant along the integral curves of the geodesic vector 
field in r(Ar\H„). In order to get a smooth first integral let 

Fkip,Y) = e-Yfz^ipX?Fk{p,Y) 

and let us define 

Fk{Kp,Y) = h{p,Y). 

So the functions F^ are smooth (non-analytic) first integrals for the geodesic flow 
on r(Ar\H„). It follows from a direct calculation that the families fzi, dAi, F 2 k-i 
and fz^ , gAi ) F' 2 k, i, k = 1,... ,n are in involution. So the geodesic flow in r(Ar\H„) 
is completely integrable in the sense of Liouville. 

Corollary 4.12. Let be the Heisenberg Lie group endowed with the standard 
metric, and let hr defined as in If Ar\H„ is the corresponding Heisenberg 

manifold, then the geodesic flow in T(Ar\H„) is completely integrable with smooth 
first integrals. 
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5. The case of a general left-invariant metric on H„ 

In this section we show how to construct first integrals for the geodesic flow on 
TH„ with arbitrary left-invariant metrics on H„. Recall that any left invariant metric 
g on is isometric to one of the form (•, ■)p defined as follows. Let (•, •) denote the 
standard inner product on ijn and let P be a symmetric positive definite operator 
on with respect to the standard inner product, which has the matrix form 



where P : D —)• 0 is symmetric and positive definite and A > 0. We can think of P 
as a symmetric matrix with positive eigenvalues. The metric (•, •)p is defined in f)„ 
by 

{X,Y)p = {PX,Y). 

Let jp{Z) be the (•, •)p-skew-symmetric operator such that 

{jp{Z)X,Y)p = {[X,YlZ)p. 

Note that 

([X,T],Z)p = {[X,Y\,PZ) = {[X,Y],XZ) 

= \{j{Z)X,Y) = X{PP-^j{Z)X,Y) 

= X{p-^j{Z)X,Y)p. 

This proves the following lemma. 

Lemma 5.1. With the assumption and notation of this section, we have that Xj{Z) = 
Pjp{Z), or equivalently, jp{Z) = XP~^j{Z), for all ^ G 3 . 

We also need the next result. 

Lemma 5.2. If f : TH„ —>• M and grad^p y) f, grad|^ y) f gradients of f with 
respect to the metrics on TH„ induced by the standard metric and {■, •)p respectively, 
then 

grad(p y) / = P grad(^ y) /, 
where P acts diagonally on T(p y)(rH„). 

Proof. By a direct calculation, 

(gradf^ y) f,S)p = (Pgradg, y) /, S) = df^^py){S) = (grad(p y) /, S) 
which proves the lemma. □ 

So, if grad(p y) / = (P, V), then 

gradfpy)/ = (p-ip,p-V). 

It follows from Proposition 1,1.,11 that fz^ is also a first integral for the geodesic flow 
on (TH,,(.,-)p). 

In order to find quadratic first integrals, we use the following lemma. 
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Lemma 5.3. An endomorphism A of V is symmetric with respect to {■,-)p if and 
only if PA = A^P, where A^ is the transpose of A with respect to the standard 
metric. 

Proof. Let X,Y E be arbitrary. Then {AX, Y)p = {X, AY)p if and only if 
{PAX, Y) = {PX, AY), and this holds if and only if {PAX, Y) = {A*PX, Y). □ 

Since P is symmetric with respect to the standard metric, there exist a basis 
Ui,... ,U 2 n of 0 such that PUi = XiUi (with A* > 0, since P is positive definite). 
This basis can be chosen orthonormal with respect to the standard metric, and 
moreover, since j{Z) acts on D as a multiple of the multiplication by on C”', 
we can assume that U 2 i = j{Zi)U 2 i-i. So we can define the operators Ai on 0 such 
that 

AiU 2 i-i = U 2 i-i AiU 2 i = U 2 i and AiUk = 0 if A: / 2z - 1, 2i. 

It follows that 

• Ai is symmetric with respect to {■, ■)p, 

• Ai commutes with jp{Z) and 

• [Ai, Aj] = 0 for all i,j. 

So Theorem 13.51 shows that we get a commuting independent family of first integrals 
for the geodesic flow on (TH^, (•, •)p). Namely 

gASP,Y) = l{AY,Y)p, 
where Ai is extended so that AiZi = 0. 

Finally, in order to obtain the remaining integrals we use the Killing fields Uf., 
k = 1,..., 2n. In fact, we have, in the same manner as in the previous section, that 

FAp, Y) = {Y, Uk)p - {jp{YAW„ Uk)p 

forms a family of first integrals for the geodesic flow such that 

{F2^-l,F2jV = 

Theorem 5.4. The geodesic flow on (TH„,(-,-)p) is completely integrable in the 
sense of Liouville. Moreover the sets 

( 1 ) 

^ = {!zA u u {F2k-i}k=i 

(2) and 

give two independent commuting families of first integrals of the geodesic flow. 

Remark 5.5. Note that the case of a general Riemannian left-invariant metric on Hn 
does not give a naturally reductive space. Also the compact quotients considered 
here are not globally homogeneous but locally. 
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